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ABSTRACT 


The  objective  of  the  present  work  was  to  experimentally 
investigate  the  effect  of  small  sinus'oidal  perturbations  on 
the  breakdown  of  circular  water  jets  issuing  from  long  tubes 
aiming  to  examine  how  the  mean  velocity  profile  for  a  fully 
developed  laminar  ,  transitional  and  turbulent  flow  as  well 
as  the  frequency  affects  the  mechanism  of  breakdown  of  the 
jet  flow. 

Observations  using  a  high  frequency  stroboscope  were 
recorded  using  a  high  speed  photographic  technique  ,  for  the 
most  distinctive  phenomena  observed,  for  three  selected  Rey¬ 
nolds  numbers  corresponding  to  laminar  ,  transitional  and 
turbulent  flow  .  Markedly  different  behavior  of  flow  break¬ 
down  and  droplet  formation  was  observed  for  each  of  these 
flows  as  frequency  waveform  or  phase  angle  was  altered. 
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I .  INTRODUCTION 


An  important  physical  phenomenon  which  has  many  practi¬ 
cal  consequences  in  technology,  ranging  from  fuel  injection 
sprays  to  ink  jet  printers,  is  the  breakdown  of  a  circular 
jet  of  fluid  issuing  into  a  quiescent  environment  . 

This  problem  has  been  extensively  studied  previously 
by  Rayleigh  and  others,  [Refs.  1,  2  and  3],  for  steady  uni¬ 
form  jets  issuing  from  an  orifice  ,  but  only  recently  have 
the  effects  of  jet  velocity  profiles  or  other  initial  condi¬ 
tions  on  the  flow  been  investigated,  [Refs.  4,  5  and  6], 

Theoretical  investigators  have  typically  assumed  jet 
breakdown  to  be  physically  similar  to  transition  from  lami¬ 
nar  to  turbulent  flow  and  have  examined  the  phenomena  using 
classical  stability  theory. 

Rayleigh  ,  [Ref.  1],  studied  the  inviscid  stability  of 
Poisseuille  flow  and  concluded  that  the  flow  was  stable  to 
infinitesimal  disturbances.  This  work  attracted  the  interest 
of  many  investigators  and  sixty  years  later  Tatsumi  [Ref.  2] 
published  his  work  related  to  the  study  of  the  velocity  dis¬ 
tribution  in  the  axisymmetric  laminar  inlet  flow  through  a 
circular  tube,  under  the  assumption  of  '  almost  similarity  ' 
of  velocity  profiles  .  He  obtained  fairly  good  agreement  of 
theoretical  with  experimental  velocity  distribution  curves 
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in  a  rather  limited  region  near  the  entrance  where  the  ap¬ 
proximation  was  valid  . 

Batchelor  and  Gill,  [Ref.  3],  have  analyzed  the  stabili¬ 
ty  of  steady  axisymmetric  jets  considering  the  typical  small 
disturbance  as  a  Fourier  component  with  sinusoidal  dependen¬ 
ce  on  both  axes  in  cylindrical  coordinates  . 

Crow  and  Champagne,  [Ref.  4]  ,  have  employed  a  loud¬ 
speaker  to  generate  a  wave  train  on  a  turbulent  jet  and 
found  that  the  phase  velocity  of  the  waves  could  be  descri¬ 
bed  by  applying  linear  theory  of  temporally  growing  insta¬ 
bility  to  a  'top  hat'  velocity  profile. 

Salwen  ,  Chester  and  Grosch,  [Ref.  5],  have  studied  the 
stability  of  Poisseuille  flow  in  a  pipe  of  circular  cross 
section  with  axisymmetric  disturbances  .  They  formulated  the 
problem  by  expanding  the  perturbation  velocity  and  pressure 
in  a  complete  set  of  orthonormal  functions  which  satisfied 
the  boundary  conditions  . 

In  a  recent  work  by  Anderson  and  Bejan,  [Ref.  6],  the 
authors  presented  a  linear  stability  analysis  of  the  large 
scale  structure  of  a  round  jet  surrounded  by  an  annular  she¬ 
ar  layer  .  Their  study  was  limited  to  the  developing  region 
near  the  jet  nozzle  at  very  large  Reynolds  number  (  Re-»-<*>). 
They  examined  the  radial  dependence  of  the  amplitude  of  gro¬ 
wing  disturbances  in  order  to  illustrate  the  extent  to  which 
the  disturbance  penetrates  into  the  jet  and  its  surroundings 


and  found  that  the  region  influenced  by  a  disturbance  was 
directly  proportional  to  the  wavelength  of  the  disturbance. 
They  also  pointed  out  that  amplified  disturbances  exhibit  a 
phase  lag  across  the  shear  layer  ,  which  may  account  for  the 
spade-like  structures  evident  in  flow  visualizations  of  tur¬ 
bulent  jets. 

In  the  present  work  we  seek  to  investigate  experimental¬ 
ly  the  effects  of  small  sinusoidal  perturbations  on  the 
breakdown  of  circular  jets  issuing  from  long  tubes  ,  which 
permit  some  control  of  their  mean  velocity  profiles  ,  for 
both  fully  developed  laminar  and  turbulent  flows. 


EXPERIMENTAL  APPARATUS  AND  PROCEDURE 


A.  DESCRIPTION  OF  EXPERIMENTAL  APPARATUS 

The  experiment  was  carried  out  in  a  vertical  falling 
jet  produced  at  the  outlet  of  a  70.0  inch  long  length  of 
0.1875  inch  inside  diameter  stainless  steel  tube  .  The  inlet 
of  the  tube  was  fed  from  a  phenum  chamber  machined  into  an 
aluminum  block  ,  which  also  contained  a  75  Watt  loudspeaker 
driver  unit  mounted  properly  to  this  block,  used  to  introdu¬ 
ce  the  pressure  perturbations  .  A  schematic  diagram  of  the 
experimental  apparatus  is  illustrated  in  Fig.  1  .  Details  of 
the  aluminum  block  and  adapter  are  shown  in  Fig.  2. 

Pressure  fluctuations  as  well  as  mean  pressure  within 
the  phenum  chamber  was  monitored  by  the  Statham  differential 
pressure  transducer  shown  in  Fig.  1  and  a  mercury  manometer. 
Mean  water  flow  was  controlled  with  a  needle  valve  and  mea¬ 
sured  with  a  0.6  gpm  rotameter  as  shown  in  Fig.  1.  In  order 
to  minimize  effects  of  geometry  at  the  exit  end  of  the  steel 
tube,  it  was  very  carefully  machined  to  ensure  sharp  edges 
and  an  absence  of  burrs. 

The  loudspeaker  driver  was  excited  with  a  Wavetec  audio 
oscillator  capable  of  producing  sine  waves  ,  as  well  as 
triangular,  unit  functions  and  ramp  functions  from  1  Hz  to 
10  MHz,  amplified  by  a  Hewlett  Packard  power  amplifier. 


1  Wavetee  H.F  VCG  Generator  ,  Model  142 

2  Hewlett  Packard  Power  Amplifier  ,  Model  467A 

3  Tektronix  Dual  Beam  Oscilloscope  ,  Type  551 

4  High  Frequency  Stroboscope  ,  Model  839 

5  Cimcon  Digital  Multimeter  DMM  51 

6  Hewlett  Packard  RMS  Voltmeter  ,  Model  3400A 

7  Intronics  Power  Supply  +15  VDC 

8  Calibration  Bridge 

9  110  VAC  Outlet 

10  Rotameter  0.6  gpm 

11  Statham  Pressure  Transducer  PM131  TC 

12  University  Sound  Driver  Unit  ID-75 

13  Mercury  Manometer 

14  Aluminum  Block 

15  Needle  Valve 

16  L.E.D  Strobolight 


Figure  1.  Schematic  Diagram  of  Experimental  Apparatus 
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Excitation  voltage  as  well  as  the  output  of  the  plenum  cham¬ 
ber  pressure  transducer  were  monitored  on  a  two  channel  Tek¬ 
tronix  oscilloscope. 

Visual  observations  of  the  jet  breakdown  process  was 
made  with  the  aid  of  a  high  frequency  LED  stroboscope,  deve¬ 
loped  at  the  Max-Planck  Institute  for  Fluids  Research1  which 
produces  a  constant  mean  light  output  from  100  Hz  to  1  MHz. 
This  device  not  only  may  be  set  for  any  phase  angle  but  also 
contains  circuitry  which  produces  automatic  phase  angle 
sweep  at  an  adjustable  rate  which  permits  periodic  phenomena 
to  be  visualized  as  they  develop.  Unfortunately  the  light 
output  was  not  sufficient  for  photographic  recording. 

Photographic  recording  of  the  observed  phenomena  was 
made  with  the  illumination  from  a  General  Radio  Microflash 
unit,  which  altough  brighter  than  the  LED  stroboscope,  lac¬ 
ked  the  sophisticated  triggering  circuitry  to  permit  control 
of  the  phase  angle  .  A  4x5  inch  view  camera  fitted  with  a 
210  mm  objective  was  used  to  photograph  on  Kodak  Tri-X  film 
which  was  force  developed  to  achieve  an  effective  film  speed 
of  ISO  3200  . 

An  overall  view  of  the  experimental  apparatus  concern¬ 
ing  and  the  photographic  equipment  used  is  shown  in  Fig.  3. 

JThe  author  extends  his  gratitude  to  Dr.  G.E.A  Meier 
who  made  the  instrument  available  for  this  work. 


17 


T 


B.  EXPERIMENTAL  PROCEDURE 

The  first  part  of  the  experiment  was  devoted  to  calibra¬ 
tion  and  determination  of  the  steady  flow  behavior  of  the 
jet  without  the  introduction  of  external  perturbations. 
As  the  flow  rate  was  increased  a  laminar  flow  column  was  ob¬ 
served  to  develope  which  converged  and  began  to  breakdown 
into  a  turbulent  flow  at  100  to  125  diameters  downstream  of 
of  the  tube  exit  corresponding  to  Reynolds  numbers  of  880  to 
2200  . 

Increasing  the  Reynolds  number  above  2200  resulted  in 
a  turbulent  jet  which  was  sudden  and  easily  identified  with 
the  naked  eye  .  Breakdown  of  the  turbulent  jet  was  observed 
with  the  high  frequency  stroboscope  and  was  observed  to  ran¬ 
ge  from  55  diameters  downstream  at  Reynolds  number  of  2288 
to  25  diameters  as  the  Reynolds  number  was  increased  to  4400 

Based  on  these  preliminary  results  ,  a  subcritical  Rey¬ 
nolds  number  of  1760  was  adopted  for  laminar  flows  and  2640 
for  turbulent  flows  and  the  effects  of  small  perturbations 
on  the  flow  at  these  two  Reynolds  numbers  was  investigated 
over  the  frequency  range  of  0  to  2000  Hz  using  the  high 
frequency  stroboscope  to  observe  the  flow. 

Based  on  these  observations  ,  photographic  records  were 
made  at  selected  frequencies  at  which  significant  characte¬ 
ristic  phenomena  occured  .  Because  of  the  difficulty  in 
accurately  describing  the  pressure  fluctuations  actually 
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produced  in  the  phenum  chamber,  oscilloscope  traces  of  these 
pressure  variations  as  well  as  the  driving  signal  were  also 
recorded.  These  results  are  presented  and  discussed  in  Chap¬ 
ter  IV.  A  discussion  of  the  effects  of  plenum  pressure  vari¬ 
ation  on  the  flow  at  the  exit  of  a  long  tube  is  presented  in 
Chapter  III  .  An  analysis  of  the  dynamic  response  of  the 
pressure  transducer  and  its  calibration  curves  may  be  found 
in  Appendices  C  and  D. 
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III.  THEORETICAL  ANALYSIS 


A.  BREAKDOWN  OF  LAMINAR  PIPE  FLOW 

One  of  the  most  important  problems  in  fluid  mechanics, 
one  which  has  attracted  the  interest  of  investigators  for 
many  years  ,  is  that  of  the  transition  from  laminar  to  tur¬ 
bulent  flow. 

The  theoretical  analysis  of  the  stability  of  laminar 
flow  is  based  upon  the  assumption  that  turbulence  is  produ¬ 
ced  by  the  growth  of  what  are  initially  small  perturbations 
or  disturbances  of  the  basic  ,  laminar  velocity  profile. 

The  theoretical  analyses  have  generally  been  limited 
to  infinitesimal  disturbances  for  which  the  differential 
equations  of  the  disturbance  can  be  linearized. 

The  basic  problem  of  the  small-disturbance  theory  is 
to  determine  whether  a  small  disturbance  is  amplified  or 
damped.  If  any  small  disturbance  is  amplified,  it  is  assumed 
that  a  transition  to  turbulent  flow  will  ultimately  occur. 
However  ,  the  small  disturbance  theory  cannot  reveal  any  of 
the  details  of  the  actual  transition  because  it  applies 
only  while  the  disturbance  is  small  .  The  question  of  how 
’small'  the  disturbance  must  be  in  order  to  apply  the  line¬ 
arized  equations  has  not  been  adequately  answered. 

An  axisymmetric  flow  is  always  stable  to  a  disturbance 
which  has  only  a  tangential  component  of  velocity,  [Ref.  7], 


Therefore  the  theoretical  analyses  have  generally  been 
limited  to  axisymmetr ical  disturbances  .  Such  analyses  have 
been  carried  out  by  many  investigators,  [Refs.  2,3  and  7]  , 
altough  some  doubt  still  remains  concerning  the  details  of 
the  mathematics. 

Steady  flow  in  a  tube  has  been  found  to  be  stable  to 
small,  axisymmetr ical  disturbances  by  a  number  of  investiga¬ 
tors, [Refs.  5,7  and  8],  and  may  be  considered  as  an  accepted 
conclus ion. 

Several  investigators  have  attempted  to  derive  a  cri¬ 
tical  Reynolds  number  ,  for  a  flow  subjected  to  external 
disturbances  ,  after  which  transition  to  turbulence  occurs. 
Leite,  [Ref.  9],  studied  the  behavior  of  small,  axisymmetri- 
cal  disturbances  in  steady  flow  in  a  tube  .  He  found  that 
they  were  always  damped  for  Reynolds  number  up  to  13000  ,  as 
long  as  their  amplitude  was  not  too  large.  Dar 1 i ng , [Ref . 10 ] , 
working  with  pulsating  flows  in  tubes  ,  found  that  the  tran¬ 
sition  Reynolds  number  dropped  from  2510  in  steady  flow  to 
1500  in  pulsating  flow.  Sarpkaya , [Ref .  11],  determined  expe¬ 
rimentally  the  critical  Reynolds  number  for  pulsating  Poi- 
seuille  flow  as  a  function  of  a  frequency  parameter  and 
velocity  ratio  .  Gilbrech  and  Combs,  [Ref.  12],  using  water 
as  working  fluid,  found  a  critical  Reynolds  number  which  was 
higher  than  2220  ,  under  some  conditions,  for  steady  flow  . 


They  also  pointed  out  that  earlier  investigators  had  found 
no  upper  limit  to  the  upper  critical  Reynolds  numbers. 

These  experimental  results  confirm  the  conclusion  that 
steady  ,  fully  developed  flow  in  tube  is  stable  to  small  di¬ 
sturbances  ,  and  the  critical  Reynolds  number  is  a  function 
of  the  frequency  of  external  disturbances  and  the  physical 
properties  of  the  fluid  .  In  spite  of  this  confirmation  , 
transition  to  turbulent  flow  is  invariably  observed  to  occur 
at  sufficiently  high  Reynolds  number  and  it  is  necessary  to 
consider  this  fact. 

The  most  obvious  explanation  is  that  transition  is 
caused  by  finite  disturbances  which  are  not  damped  ,  even 
though  smaller  disturbances  are  .  This  possibility  is  stron¬ 
gly  supported  by  the  observation  made  by  Gilbrech  and  Combs 
[Ref,  13],  that  the  critical  Reynolds  number  decreases  as 
the  magnitude  of  the  disturbance  increases  . 

The  theoretical  analysis  of  finite  disturbances  has 
been  limited  by  mathematical  difficulties  ,  however  ,  Stuart 
[Ref.  14],  has  shown  that  a  flow  may  be  unstable  to  finite 
disturbances  ,  while  it  is  stable  to  infinitesimal  ones.  In 
addition,  Spielberg  ,  [Ref.  15],  has  pointed  out  that  stabi¬ 
lity  to  two  dimensional  disturbances  does  not  necessarily 
imply  stability  to  three  dimensional  disturbances.  This  lat¬ 
ter  result  creates  a  problem  in  comparing  accurately  the 
theoretical  and  experimental  results  . 


In  the  present  work,  it  was  assumed  that  an  axisymmetri- 
cal  disturbance  was  propagated  without  confirming  the 
non-existance  of  three  dimensional  disturbances  . 

Another  explanation  of  the  origin  of  turbulence  in  pipe 
flow  may  be  that  the  transition  occurs  in  the  entrance 
region  of  the  pipe,  in  which  the  velocity  distribution  is  of 
the  boundary-layer  type  .  A  theoretical  analysis  by  Tatsumi, 
[Ref.  2],  indicates  that  the  flow  in  the  entrance  region  is 
unstable  to  small  disturbances.  He  found  a  point  of  maximum 
instability  near  the  entrance  to  the  tube,  where  the  maximum 
Reynolds  number  for  stability  was  19400  .  However  ,  the  cri¬ 
tical  Reynolds  number  increases  rapidly  in  both  directions, 
so  that  one  cannot  conclude  that  this  local  instability  will 
lead  to  transition.  This  is  because  amplification  required 
to  produce  transition  remains  undetermined. 

Sparrow  and  Lin,  [Ref.  16],  pointed  out,  in  their  effort 
to  determine  the  detailed  nature  of  the  flow  development  in 
the  entrance  region,  that  even  for  a  laminar  flow, the  velo¬ 
city  problem  for  the  entrance  region  does  not  yield  an  exact 
solution,  regardless  of  the  shape  of  the  duct  cross  section. 
Their  explanation  for  the  difficulties  in  the  analysis  was 
that  the  essential  nonlinearity  of  the  inertia  terms  which 
appear  in  the  equation  of  motion,  preclude  linear  analysis. 

From  a  more  sophisticated  point  of  view  ,  the  problem 
becomes  especially  complicated  if  we  attempt  to  reveal  the 


nature  of  motion  in  the  non-linear  range  of  boundary  layer 
instability  and  the  onset  of  turbulence. 

Klebanoff  and  others  ,  [Ref.  17]  ,  who  attacked  the  pro¬ 
blem  in  the  general  form  ,  pointed  out  that  despite  the  suc¬ 
cess  of  the  linearized  theories  in  revealing  the  nature  of 
the  initial  stages  of  boundary-layer  instability  ,  there  re¬ 
mains  a  deep  void  in  the  undestanding  of  the  subsequent  non¬ 
linear  behavior  and  the  actual  breakdown  of  the  laminar 
boundary-layer . 

Miller  and  Fejer,  [Ref.  18]  ,  who  studied  the  transition 
in  Blasius-type  boundary-layer  produced  by  a  free  stream  ha¬ 
ving  an  oscillatory  component  of  velocity  ,  found  that  the 
transition  Reynolds  number  depends  only  on  the  amplitude  of 
of  the  oscillations  ,  and  that  the  dimensionless  transition 
length  is  a  function  only  of  the  frequency. 

Obremski  and  Fejer  ,  [Ref.  19]  ,  continued  the  previous 
work  and  pointed  out  that  below  a  critical  value  of  the  pa¬ 
rameter  [  (Re)Ng=  LAU/2rrv  ]  ,  transition  occurs  at  a  relati¬ 
vely  constant  Reynolds  number  ,  which  appeared  independent 
of  the  amplitude  and  frequency  of  the  oscillation  ,  at  least 
over  the  range  investigated. 

Finally  Landahl  ,  [Ref.  25],  used  kinematic  wave  theory 
to  determine  under  what  conditions  breakdown  of  a  steady  or 
unsteady  laminar  flow  into  high  frequency  oscillations 
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should  occur  .  Isolating  the  three  different  ingredients'  of 
the  mechanism  of  breakdown  ,  he  pointed  out: 

a.  The  primary  instability  process  itself  is  only 
incidental  to  breakdown  and  transition  ,  since  one  can 
easily  produce  a  breakdown  condition  by  a  sufficient¬ 
ly  intense  local  disturbance  of  the  shear  flow  ,  for 
example  through  free  stream  disturbances. 

b.  If  the  breakdown  is  to  be  self  maintained  as  it 
travels  downstream  however  ,  hydrodynamic  instability 
is  required.  In  a  sense  ,  the  breakdown  mechanism  is 
the  one  most  essential  in  transition  of  a  boundary- 
layer  to  turbulence  rather  than  the  classical  hydro- 
dynamic  instability  one,  since  the  former  represents 

a  strongly  irreversible  process. 

c.  An  unstable  small  •  amplitude  wave  packet  of  the 
Tolmien-Schl icht ing  type  may  amplify  for  a  while  ,  as 
it  moves  downstream  ,  but  once  it  has  passed  through 
some  streamwise  position  in  the  boundary-layer  it 
will  leave  the  shear  flow  practically  as  undisturbed 
as  before  the  passage  of  the  packet  ,  except  for  a 
small  deformation  of  the  mean  flow  of  second  order 
in  disturbance  amplitude  .  At  breakdown  on  the  cont¬ 
rary  there  is  an  irreversible  redistribution  of  the 
basic  shear  flow  vorticity  due  to  a  nonlinear  recti¬ 
fication  mechanism  . 

In  view  of  the  difficulties  encountered  in  applying 
the  small  disturbance  theory  to  steady  axisymmetr ical  flow  , 
it  is  likely  that  considerable  time  and  effort  will  be  re¬ 
quired  to  obtain  an  exact  solution  for  unsteady  flow  pro¬ 
blems  ,  especially  if  we  consider  a  three  dimensional  flow  . 

The  use  of  a  ’  quasi-steady  ’  analysis  would  reduce  the 
problem  to  one  which  could  be  handled  more  easily.  Using 
this  technique  one  would  analyze  an  instantaneous  velocity 
profile  as  though  it  were  a  steady-flow  problem. 


By  considering  the  stability  of  the  velocity  profiles 
at  various  times,  one  might  then  draw  some  conclusions  about 
the  stability  of  the  unsteady  flow  an  approach  used  by  Shen, 
[Ref.  21],  as  well  as  Greenspan  and  Benney,  [Ref.  22]. 

In  the  present  work,  although  Reynolds  Numbers  as  high 
as  19400  were  not  reached  ,  we  recognize  that  the  flow  in 
the  entrance  region  has  ,  in  all  probability  ,  affected  the 
downstream  flow  and  may  well  have  introduced  turbulence  sub¬ 
sequently  damped  along  the  tube  ,  for  low  Reynolds  numbers  , 
creating  a  laminar  or  laminar-like  flow  at  the  tube  exit. 


B.  NONSTEADY  FLOW  IN  A  CIRCULAR  TUBE 

In  order  to  study  how  the  pressure  disturbances,  intro¬ 
duced  at  the  entrance  of  a  long  tube  used  to  form  the  jet  , 
affect  the  flow  at  the  exit  of  the  tube  ,  which  constitutes 
the  initial  conditions  of  the  flow  of  interest  ,  we  may  exa¬ 
mine  the  nonsteady  flow  in  a  long  tube  of  circular  cross 
section  .  We  adopt,  as  shown  in  Fig.  4  ,  cylindrical  coordi¬ 
nates  whose  x-axis  is  identified  with  the  center  line  of  the 
pipe. 

The  equation  of  continuity  is: 


V(pV)  +  —  =  0 
*t 


(Eq.  1) 


In  cylindrical  coordinates  Eq.  1  becomes: 
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(Eq.  2) 


1  ^  1  0 

7 rP  <rf>Vr)  +7^<pV9>  +  >i <pVx>  +>«  ’ 


x - 


Figure  4.  Cylindrical  Coordinates  in  a  Round  Pipe. 


Assuming  incompressibility,  (o  =  constant),  and  a  non¬ 
twisting  flow  (vQ  =  0),  Eq.  2  may  be  written  as: 

1  ^  ^ 

<rPvr>  +  r-  <Pvx>  =  0 

pop  Ox 


^v_  V 

J  - 5  +  — l  +  -I  =  0 
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(Eq.  3) 


Neglecting  external  forces  ,  (i.e  gravity)  ,  the  Navier- 
Stokes  equations  in  two  components  can  be  written: 
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(Eq.  5) 


If  we  assume  quasi-steady  flow,  the  equation  of  conti¬ 
nuity  Eq.  3  becomes: 


^  v 

- -  =0  or :  v  =  v( r  ,  t ) 

dt 


(Eq.  6) 


This  indicates  that  the  velocity  in  the  direction  pa¬ 
rallel  to  the  center  line  is  a  constant  at  each  moment. 
Substituting  =  0  into  Eq.  5  we  get: 


dP 


=  0  or:  P  =  P(x,t) 


(Eq.  7) 


Inserting  these  results  in  the  equation  of  motion 
(Eq.  4)  : 
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(Eq.  8) 


where:  V  =  — 


i s  the  kinematic  viscocity  of  the  fluid. 


29 


C.  SOLUTION  FOR  A  UNIFORM  PULSATING  FLOW 


Assuming  that  the  pressure  pulses  introduced  by  the 
drive  unit,  (a  loudspeaker  driver),  are  propagated  instantly 
along  the  flow  direction  ,  (  justified  by  the  fact  that  the 
mean  flow  velocity  is  very  small  compared  with  the  speed  of 
sound  in  the  fluid  )  ,  with  the  aid  of  Eqs.  6,  7  and  8  we 
can  show  that  the  pressure  gradient  ,  (t>P/3x  ),  becomes  a 
function  only  of  time  ,  t. 

Using  Fourier  series  we  can  express  the  pressure  gra¬ 
dient  as  follows  : 


1  ^P  00  00 

-T-  =  “o  +  iS  b„sinfit) 

pox 


(Eq.  9) 


or,  for  the  sake  of  simplicity  of  calculation  ,  in  a  complex 
form: 


1  hP  ~ 

=  %  +  n=l*nein  •  tfn  =  °S  ”  ibn 


(Eq.  10) 


where  :  ,  bR  are  constants  representing  the  amplitudes  of 

elemental  vibrations. 

The  corresponding  solution  for  the  axial  velocity  may 
be  assumed  to  be  of  the  form  : 


vx  =  vx0 


+  vnicos(pt)  +  jffi  Vp2  s  i  n(nt) 


(Eq.  11) 
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or,  in  complex  form: 


vx  *  vx  +  ^  vneini  .  vn  =  vn  -  ivn 


(Eq.  12) 


where  coefficients  v  and  v  are  assumed  to  be  functions 

ni  n* 

of  the  radial  coordinate  ,  r  ,  only. 

Substituting  Eqs.  10  and  12  into  Eq.  8  and  comparing 
the  terms  of  the  same  order  ,  we  obtain  the  following  Dif¬ 
ferential  Equations  for  the  coefficients  as  independant 
variables. 


dr2  r  dr  v 


(Eq.  13) 


and : 
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dr2  r  dr  v  n  v 


(Eq.  14) 


After  integration  we  get: 
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(Eq.  15) 


where : 


and : 


A„  , B„  =  constant 
o  o 
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(Eq.  16) 


where:  5  =  V  (n/v) 

and  Jo  and  Ko  are  the  Bessel  functions  of  the  first  and 
second  kind  of  scro  order  respectively. 

The  solution  for  the  x-wise  component  of  velocity  is 
given  by  the  real  part  of  the  following  expression: 


v  - - -  A0lnr  +  B0  + 

v  4 
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On  the  axis  of  the  pipe  ,  (  r  =  0  )  ,  ln(r)  =  and 

Ko(Sri1/2)  =  oo  ,  since  vx  is  finite  ,  we  conclude  that: 

A  =  E  =  0  (£<1*  18^ 

o  n 

The  Boundary  Condition  at  the  wall  is  given  by: 
vx  =  0  at  r  =  R  (Eq.  19) 

where:  R  is  the  radius  of  the  pipe. 

Inserting  this  boundary  condition  ,  Eq.  19  ,  into 
Eq.  17  we  get: 


B 


o 


ill* 

n=1  lin 


+  DnJo 


(SRi3/2> 


(Eq.  20) 


Inserting  ,  Eqs.  18  and  20  into  Eq.  17,  we  finally  get 
the  velocity  component  vx  l 
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(Eq.  21) 
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The  total  mean  mass  flow  ,  m  ,  which  is  transmitted  in 
the  x-direction  is  given  by  the  integral  over  time  and  ra¬ 
dius  of  the  axial  velocity  component: 


l  2  rr  R 

m  =  - J  J2  Tlvxrcjr  = 
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(Eq.  22) 


Inserting  the  time  mean  pressure  gradient: 

f  Til 
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(Eq.  23) 


we  get: 


ttR4 

t>P 

8p 

(Eq.  24) 


It  is  seen  that  the  total  mean  mass  flow  ,  in  pulsating 
motion  ,  is  identical  to  that  of  a  steady  Poisseuille.  flow 
with  the  same  pressure  gradient,  as  the  mean  pressure  gra¬ 
dient  in  the  pulsating  flow. 

A  mean  velocity  ,  U  ,  may  be  defined  by: 

m 

u  =  - ;  (Eq.  25) 


Combining  Eqs.  24  and  25  we  get: 


R2  I  dP 
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(Eq.  26) 


We  may  express  the  x-wise  velocity  as  the  sum  of  a 

perturbation  velocity,  f  and  the  mean  stream  velocity  of 

the  fluid.vw„  • 

*  xs 


vx  =  vxs  +  vx 


(Eq.  27) 


Introducing  the  Hagen-Poiseui 1 le  relation  for  v  We 

X  s 


get : 


=  2  1- 


(Eq.  28) 


Dividing  Eq.  27  by  the  total  mean  velocity  ,  U  ,  we 
get  the  non-dimensional  expression: 
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(Eq.  29) 


Combining  Eqs.  21  and  29  we  get  for  the  perturbation 


velocity  ,  vx  • 
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where: 


A  = 


ber (SR)ber (Sr )  +  be i (SR)be i (Sr ) 


B  = 


ber ^ (SR)  +  be i ^ (Sr) 
bei (SR)ber (Sr )  -  ber (SR)bei (Sr ) 


ber^ (SR)  +  be i ^ (SR) 


(Eq.  31) 


where:  ber  and  bei  are  the  Kelvin  functions  . 

With  the  aid  of  Eqs.  10  and  26  we  can  obtain  the  cor¬ 
responding  pressure  gradient  relation: 
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(Eq.  32) 


where : 

Re  =  2RU  /  v 

Introducing  the  friction  factor  for  laminar  flow  : 

f  =  64  /  Re 

into  Eq.  32,  we  get  : 

J)P  pO^  r  tt’  ot  tt1  bn  i 

"  ~  =  ~Z~f  1  *  ith — cosfot)  +  — sinfit)  (Eq.  33) 
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*  n/a0  and  bn/ofo  are  dimensionless  ratios  of 
amplitude  of  periodic  variation  of  pressure  gradient  norma¬ 
lized  with  the  mean  amplitude  . 


D.  VELOCITY  DISTRIBUTION  AND  ASYMPTOTIC  EXPRESSIONS 


The  velocity  profiles  for  laminar  ,  pulsating  flow  cons¬ 
titute  the  initial  conditions  for  any  theoretical  analysis 
of  stability.  Although  no  exact  theoretical  analysis  of  sta¬ 
bility  is  presnted  in  this  work,  solution  for  the  velocity 
profiles  will  serve  to  identify  the  dimensionless  parameters 
which  are  significant  and  will  permit  some  qualitative  de¬ 
duction  concerning  the  stability  of  pulsating  flow. 

From  Eq.  16  we  found  that  the  periodic  part  of  the  velo¬ 
city  profile  is  characterized  by  the  parameter  Sr  = /n/v  R. 

The  distribution  of  a,,/  ocQ  ,  the  dimensionless  velocity 
amplitude,  obtained  by  Gilbrech  and  Combs  ,  [Ref.  16]  ,  is 
plotted  in  Figure  5. 

Uchida,  [Ref. 23],  obtained  solutions  for  limiting  values 
of  the  parameter  oR  by  means  of  asymptotic  expansions  of  the 
Bessel  function: 


1 .  Slow  Oscillation  (  SR  <<  1  ) 

When  a  highly  viscous  fluid  oscillates  slowly  in  a 
narrow  pipe  ,  oR  becomes  small  and  we  can  approximate: 


ber(SR)  —  1  and  bei(SR)  —  0 

Using  this  approximation,  we  get  from  Eqs.  28  and  29 
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Figure  5.  Oscillation  Amplitude  of  the  Local  Velocity 
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(Eq.  34) 


From  this  we  conclude  that  the  velocity  distribution 
is  a  parabola,  as  in  the  case  of  steady  Hagen  -  Poisseuille 
flow,  while  the  magnitude  of  velocity  varies  periodically  in 
phase  with  that  of  the  pressure  gradient. 

2.  Rapid  Oscillation  (  5r  -»-<*>  ) 


When  a  fluid  of  low  viscocity  oscillates  rapidly  in 
a  large  pipe  ,  then  parameter  oR  becomes  large. 

Assuming  §R  >  10  ,  asymptotic  expansion  of  Bessel 
functions  may  be  introduced: 


eSR/2  r  r£R  n  r  1 
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eSR/2  f  rSR  n  1  1 

be  i  (SR)  =  - lsin - +  0 

2  tt5R  12  8  OR 


(Eq.  35) 


To  examine  the  motion  of  the  fluid  near  the  center 


of  the  pipe  ,  we  put: 


Sr  oo  and  5r  — *■  0 


As  before  we  get: 
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(Eq.  36) 

From  this  result  we  may  conclude  that  in  the  case  of 
rapid  oscillation  ,  fluid  flows  at  the  center  of  pipe  with 
the  phase  lag  of  90°  behind  the  wave  of  pressure  gradient 
and  its  amplitude  diminishes  with  increasing  frequency. 

To  examine  motion  of  the  fluid  near  the  wall  of  the 
pipe  ,  we  put: 

5R  00  and  8r  —•*«*> 

Using  the  same  procedure  we  get: 


— =2|1-— 
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[AfA  •  {sinM.£ 

l  R2J  n  «  ( SR ) 2  l  » r 


S(R-r) 


'  r  5  ll 
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ocQ  ( 3TR ) 2 


f  P" 

<  -cos(n  t)+  w— e 


5(R-r ) 


'cos  nt - (R-r)  >  (Eq.  37) 

i/2  $  JJ 


Because  of  the  obvious  complexity  of  the  expression, 
it  is  difficult  to  draw  any  simple  conclusions  other  than  to 
note  that  a  local  maximum  in  velocity  developes  in  the  neigh 
borhood  of  the  wall  as  it  is  shown  in  Fig.  5  .  We  can  also 
get  a  better  confirmation  from  Figs.  6  ,  7  and  8  ,  obtained 
from  [Ref.  23],  for  values  of  5r  =  1 ,  3  and  10  respectively. 
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E.  SECTIONAL  MEAN  VELOCITY 


Since  the  instantaneous  mass  flow  and  accordingly  sec¬ 
tional  mean  velocity  ,  are  periodic  functions  ,  we  seek  to 
link  the  mean  velocity  with  the  applied  pressure  gradient* 
Sectional  mean  velocity  ,  denoted  by  v^^,  is  given  by: 

R 


xm 


1  r 


vxrdr 


(Eq.  38) 


Introducing  Eq.  34  and  integrating  we  get: 


xm  =  i  +  ^  fs  8 
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where: 


(Eq.  39) 


dfber(SR)]  ,  d(ber(SR)] 

ber'(SR)  =  - - _ -  ,  bei(SR)  =  -- - -  (Eq.  40) 


d  (Sr) 


d  (Sr) 


and: 


C  = 


D  = 


ber  (SR) be i*(SR)  -  be i  (SRjber^SR) 
be r 2 (SR)  +  be i 2 (Sr) 
ber  (Sr) be r'( Sr)  +  be  i  (SR)be  i^SR) 


(Eq.41 ) 


ber2(5R)  +  bei2(5R) 

An  example  of  simple  periodic  pulsation  is  given  by: 
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V  _\ 


represents  the  coefficient  of  phase  lag  in  the  wave  of 
pressure  gradient. 

In  Fig.  9  ,  obtained  from  [Ref.  23]  ,  are  plotted  the 
coefficients  of  amplitude  and  phase  lag  of  mean  velocity  as 
a  function  of  the  parameter  oR. 


Figure  9.  Coefficients  of  Amplitude  and  Phase 
Lag  of  Mean  Velocity. 


F.  SURFACE  FRICTION  AND  BALANCE  OF  FORCE 

The  instantaneous  frictional  force  acting  on  the  pipe 
wall  is  given  bys 


r  = 


r=R 


(Eq.  44) 


Dividing  Eq.  44  by  (  pU2/2  )  ,  and  combining  with  the 
previously  obtained  Eqs.  32  ,  39  and  40  we  can  get: 
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(Eq.  45) 


Considering  the  simple  periodic  pulsation  described  by 
Eq.  42  we  get: 
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pU2  Re  l  a0‘  Sr  Sr  J 
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where: 
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(Eq.  46) 


represents  the  coefficient  of  amplitude  of  shearing  stress 
at  the  wall  ,  and: 


*T=  Un  [7 


represents  the  coefficient  of  phase  lag  of  shearing  stress 
with  respect  to  the  wave  of  pressure  gradient. 
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In  the  Fig.  10  below  ,  obtained  from  [Ref.  23],  it  is 
shown  how  the  coefficient  of  amplitude  and  phase  lag  of 
shear  stress  varies  with  Sr  . 


Figure  10.  Coefficients  of  Amplitude  and  Phase  Lag 
of  the  Shearing  Stress 

If  we  compare  Figs.  9  and  10  we  observe  no  substantial 
differences  in  amplitude  or  phase  lags.  Moreover,  the  phase 
of  the  sectional  mean  velocity  is  greatly  delayed  with  res¬ 
pect  to  the  pulsating  wave  of  pressure  gradient,  while  that 
of  shearing  stress  is  less  delayed  from  it. 

Using  Eq.  8  and  integrating  over  the  section  of  the 
pipe,  from  r  =  0  to  r=R  ,  we  can  obtain  a  relationship  between 
the  force  and  acceleration  of  amass  of  fluid  enclosed  in  the 
circular  cylinder  of  unit  length: 


dv 

-  ttR2 — =  ttR2o — — 


+  2  ttRt 


or : 


1  -  dP/dx 

2  pU2/4R 


dYim/dt  4T 

02/2R  pO2 


(Eq.  47) 


For  the  present  example  of  a  simple  periodic  pulsation 
we  find  the  mean  frictional  force  ,  involved  in  the  shear 
stress  ,  is  just  balanced  by  the  term  of  mean  pressure  gra¬ 
dient  and  there  results  no  steady  acceleration.  The  various 
terms  appearing  in  Eq.  47  are  shown  in  Fig.  11  below  ,  ob¬ 
tained  from  [Ref.  23]. 


Figure  11.  Periodic  Force  Component  (A  -  Pressure  Gradient, 
B  -  Acceleration  ,  C  -  Shearing  Force) 


IV.  RESULTS  AND  CONCLUSIONS 


A.  PRELIMINARY  RESULTS 

After  a  preliminary  survey  of  the  phenomenon  of  jet 
breakdown  at  a  single  laminar  and  single  turbulent  Reynolds 
number,  as  mentioned  above,  it  was  found  that  the  field  was 
so  rich  in  information  that  a  third  ,  transitional  Reynolds 
number  of  2200  was  added  to  the  survey  .  Observations  repor¬ 
ted  are  of  the  most  dinstinctive  phenomena  at  each  of  the 
three  Reynolds  numbers  and  were  observed  to  occur  most  fre- 
quentlyat  frequencies  which  were  multipliles  of  33  and  50  Rz 
Results  include  two  photographs  at  each  data  point  of  two 
separate  phases  along  with  an  oscillogram  showing  both  the 
excitation  waveform  and  the  pressure  perturbation  it  produ¬ 
ced,  shown  in  Figs.  13  through  36.  Table  I  is  a  summary  of 
the  observed  breakdown  length  as  a  function  of  the  frequency 
and  Table  II  is  a  summary  of  oscilloscope  and  flow  parame¬ 
ters,  corresponding  to  each  figure. 

As  Hoyt  and  Taylor,  [Ref.  24],  pointed  out,  jets  are  un¬ 
stable  to  two  kinds  of  disturbances,  axial  and  helical.  The¬ 
ory  suggests  that  the  axial  instability  (  corresponding  to 
the  initial  waves  on  the  jet  surface  )  die  out  after  a  few 
nozzle  diameters  downstream  and  the  helical  instabilities, 
amplified  by  aerodynamic  resistance  ,  lead  to  the  final  jet 


breakup.  The  initial  slight  helical  motion  is  amplified  un¬ 
til  the  entire  jet  may  be  described  as  having  a  corksrew- 
like  motion  .  This  is  well  confirmed  by  the  photos  obtained 
in  the  present  work  of  the  jet  formation  before  and  after 
the  breakdown  and  is  in  good  agreement  with  the  observations 
obtained  by  Crow  and  Champagne,  [Ref.  4]. 


TABLE  I 


SUMMARY  OF  BREAKDOWN  DISTANCE 


(Hz) 

#d 

(Re=l76Q) 

#d 

(Re  =  2200 ) 

#d 

(Re=2640 

0 

118.0 

125.0 

30.0 

33 

112.6 

120.6 

31.5 

50 

110.0 

112.0 

34.0 

66 

106.3 

103.7 

32.2 

100 

102.0 

98.5 

29.5 

133 

98.6 

99.2 

33.7 

150 

95.5 

100.0 

35.0 

166 

92.7 

101.7 

37.8 

200 

88.0 

102.5 

41.0 

233 

85.2 

101.2 

45.2 

250 

83.5 

99.0 

49.0 

266 

81.5 

101.5 

47.6 

300 

91.2 

88.0 

45.0 

333 

96.1 

97.5 

44.1 

350 

101.8 

110.3 

43.3 

366 

110.0 

121.0 

42.0 

400 

112.0 

125.0 

42.5 

Re  =  Reynolds  Number 
f  =  Introduced  Frequency 
#d  =  Number  of  Diameters 
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TABLE  II 


SUMMARY  OF  EXPERIMENTAL  CONDITIONS 


Fig. 

#  Re 

f 

(Hz) 

vp 

(mV/cm) 

Vg 

(mV/ cm) 

t/cm 

(msec/cm 

13 

1760 

50 

5 

10 

5 

14 

1760 

100 

5 

10 

5 

15 

1760 

150 

5 

10 

2 

16 

1760 

200 

5 

10 

2 

17 

1760 

250 

5 

10 

2 

18 

1760 

270 

5 

10 

2 

19 

1760 

300 

5 

10 

2 

20 

1760' 

380 

5 

10 

2 

21 

2200 

50 

5 

10 

5 

22 

2200 

100 

5 

10 

5 

23 

2200 

150 

5 

10 

2 

24 

2200 

200 

5 

10 

2 

25 

2200 

233 

5 

10 

2 

26 

2200 

250 

5 

10 

2 

27 

2200 

266 

5 

10 

2 

28 

2200 

300 

5 

10 

2 

29 

2200 

380 

10 

10 

1 

30 

2200 

787 

20 

20 

0.5 

31 

2640 

50 

5 

10 

5 

32 

2640 

100 

5 

10 

2 

33 

2640 

200 

5 

10 

2 

34 

2640 

250 

10 

10 

2 

35 

2640 

300 

10 

10 

1 

36 

2640 

600 

10 

10 

0.5 

Re  : 

Reynolds  number 

f  s 

Introduced  Frequency 

vp  = 

Signal 

Ampl i tude 

from  Pressure  Transducer 

Vg  * 

Signal 

Amp  I i tude 

from  Frequency 

Generator 

1  I  I  I  I  I - 1 - 1 - 1 


0.0  50.0  100.0  150.0  200.0  250.0  300.0  350.0  400 

FREQUENCY  (HZ) 

Figure  12.  Breakdown  Length  Versus  Frequency 
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1.  Laminar  Flow  at  Reynolds  Number  1760 

At  this  Reynolds  number  the  unexcited  flow  started 
to  breakdown  at  approximately  118  diameters  from  the  tube 
exit.  Increasing  the  frequency  to  266  Hz  produced  a  gradual¬ 
ly  movement  of  the  breakdown  point  upstream  to  approximately 
81.5  diameters  from  the  exit.  Above  this  frequency  a  rever¬ 
sed  trend  appeared  and  the  breakdown  point  started  to  move 
downstream  as  the  flow  became  laminar  similar  to  the  unexci¬ 
ted  flow  ,  the  only  difference  being  in  the  droplets  formed 
after  the  breakdown.  This  phenomenon  may  be  characterized  as 
a  re  laminar i zat ion  of  the  flow  and  was  almost  complete  at 
380  Hz.  Further  increasing  of  the  frequency  up  to  1000  Hz 
produced  no  significant  effect  in  the  jet  flow. 

Results  for  Reynolds  number  1760  are  shown  in  Figs. 
13  through  20  .  As  can  be  observed  from  the  photos  ,  break¬ 
down  point  is  not  fixed  and  oscillates  somewhat  upstream  and 
downstream  ,  creating  a  breakdown  region  .  The  amplitude  of 
this  oscillation  was  reduced  as  the  perturbation  frequency 
was  increased  ,  starting  from  about  14  diameters  and  approa- 
approching  4  diameters  as  the  frequency  changed  from  200  to 
266  Hz  .  At  266  Hz  an  especially  interesting  phenomenon  was 
observed  concerning  the  breakdown  point  and  the  droplet  for¬ 
mation  .  The  breakdown  of  the  jet  occured  at  the  minimum  ob¬ 
served  distance  from  the  tube  exit  at  81.5  diameters. 


..... 


The  oscillation  of  the  breakdown  point  was  minimized  to  a- 
bout  4  diameters  and  the  droplets  formed  ,  as  was  observed 
with  the  aid  of  the  stroboscope  ,  were  exactly  the  same  and 
appeared  stationary  in  space.  This  is  shown  in  the  jet  pho¬ 
tos  of  Fig.  18. 

A  plot  of  the  breakdown  distance  as  a  function  of 
frequency  is  provided  in  Fig.  13  where  we  see  that,  in  the 
frequency  range  from  zero  to  266  Hz,  an  almost  linear  varia¬ 
tion  of  the  breakdown  point  with  frequency  occurs.  The  mini¬ 
mum  at  266  Hz  may  be  characterized  as  a  critical  frequency 
for  this  specific  Reynolds  number  .  Above  this  frequency  one 
also  observes  an  almost  linear  variation  of  the  average 
breakdown  point  to  380  Hz.  This  region  may  be  characterized 
as  relaminarization  region.  The  slope  of  the  breakdown  curve 
is  greater  than  at  lower  frequencies  although  increasing  the 
frequency  above  380  Hz  results  a  zero  slope  curve  ,  as 
discussed  above. 

2 .  Transitional  Flow  at  Reynolds  Number  2200 

At  this  Reynolds  number  the  unexcited  flow  started 
to  breakdown  at  approximately  125  diameters  from  the  tube 
exit  .  Increasing  the  frequency  to  100  Hz  ,  the  breakdown 
point  moved  upstream  up  to  about  98  diameters  from  the 
tube  exit  .  The  oscillation  of  the  breakdown  point  was 
reduced  nt  this  case,  from  12  to  6  diameters. 
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Increasing  the  frequency  to  266  Hz  resulted  in 
only  slight  variation  of  the  breakdown  point  .  The  only 
unusual  observation  in  this  range  of  frequency  was  the 
drastic  reduction  of  the  amplitude  of  breakdown  point 
oscillation  as  the  frequency  was  increased  . 

An  interesting  pressure  effect  was  observed  at  266 
Hz,  as  can  be  seen  at  Fig.  27  ,  where  a  maximum  amplitude 
of  pressure  signal  was  recorded. 

An  interesting  change  was  observed  when  the  fre¬ 
quency  reached  300  Hz  .  The  breakdown  point  was  displaced 
suddenly  about  25  diameters  upstream  and  nearly  stabilized 
at  88  diameters  .  The  droplets  viewed  with  the  aid  of  the 
stroboscope  appeared  immobile  and  any  oscillation  of  the 
breakdown  point  was  difficult  to  observe  .  Based  on  a  series 
of  high  speed  photos  it  was  concluded  that  the  breakdown 
point  was  oscillating  with  an  amplitude  of  about  2  diameters 
Also  observed  at  this  frequency  that  the  waveforms  obtained 
from  the  signal  generator  and  the  pressure  transducer  were 
almost  in  phase  and  moreover,  the  latter  waveform  appeared 
nearly  sinusoidal  as  shown  in  Fig.  28. 

Further  frequency  increases  resulted  in  a  flow  rela- 
minarization  as  was  observed  and  discussed  in  the  previous 
case.  A  completely  laminar  flow  was  observed  at  a  frequency 
of  380  Hz  and  no  additional  significant  effects  were  obser¬ 
ved  up  to  2000  Hz  . 


63 


Figure  21.  Photos  Illustrating 
Trace  of  Pressure  a 


Figure  28. 


Photos  Illustrating  Jet  Breakdown  and  Oscilloscope 
Trace  of  Pressure  at  Be»J200  and  f=30Q  Hz 
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Refermg  to  Fig.  12  ,  we  may  draw  some  conclusions 
concerning  the  average  breakdown  point.  Up  to  100  Hz  a  gra¬ 
dually  movement  of  the  breakdown  point  upstream  is  observed 
which  is  nevertheless  more  rapid  than  the  laminar  case  as 


can  be  seen  from  comparison  of  the  slope  of  the  curves 
within  this  region  .  A  nearly  stable  region  then  appears  in 
the  frequency  range  between  100  Hz  and  266  Hz  where  the 
breakdown  point  is  nearly  fixed  at  100  diameters.  This  re¬ 
gion  may  be  characterized  as  a  stabilized  region  and  was 
not  observed  at  any  other  Reynolds  number.  A  singular  point 
was  also  observed  at  300  Hz  as  discussed  previously.  This 
frequency  may  be  characterized  as  a  critical  frequency  for 
the  specific  Reynolds  number  .  Finally  ,  we  can  see  from 
Fig.  12  that  the  relaminar izat ion  rate  was  also  greater  for 
transitional  flow,  producing  a  steeper  slope  of  the  curve  in 
the  frequency  range  from  300  Hz  to  400  Hz. 


3.  Turbulent  Flow  at  Reynolds  Number  2640 


At  this  Reynolds  number  it  was  observed  that  the 
unexcited  flow  began  to  breakdown  at  an  average  distance  of 
30  diameters  and  no  significant  effect  on  the  jet  was  obser- 
rved  up  to  a  frequency  of  100  Hz.  Increasing  the  frequency  a 
gradual  movement  of  the  average  breakdown  point  downstream 
observed,  which  reached  a  maximum  of  50  diameters  at  250  Hz. 


Further  increase  of  the  frequency  produced  a  very 
slow  movement  of  the  average  breakdown  point  upstream  to 
about  42  diameters  at  a  frequency  of  366  Hz.  The  breakdown 
point  then  remained  almost  unchanged  up  to  2000  Hz  .  The 
only  further  change  observed  in  the  high  frequency  range  , 
from  366  Hz  to  2000  Hz,  was  the  appearance  of  a  corksrew- 
like  motion  of  the  jet  which  became  more  pronounced  with 
increasing  frequency  ,  making  more  readily  observable  the 
helical  mode  ,  which  dominates  the  breakdown  mechanism  . 
These  results  are  shown  in  Figs.  31  through  36. 


B.  CONCLUSIONS. 

From  the  results  ,  discussed  above  ,  we  may  draw  the 
following  conclusions  concerning  the  effect  of  small  pertur¬ 
bations  on  jet  breakdown: 

a.  In  the  case  of  laminar  flow,  (Re=1640)  ,  the  breakdown 
point  moves  upstream  as  the  frequency  is  increased  up 
to  266  Hz.  Downstream  of  the  breakdown  point  the  jet 
flow  becomes  turbulent  and  breaks  into  discrete  dro¬ 
plets.  Above  266  Hz  this  effect  is  observed  to  reverse 
itself  and  the  jet  returns  to  its  laminar  nature. 

b.  In  the  case  of  transitional  flow  ,  (Re=2200)  ,  the 
breakdown  point  moves  initially  upstream  as  frequency 
is  increased  to  100  Hz  .  In  the  range  of  100  Hz  to 
266  Hz  an  almost  stable  behavior  of  the  breakdown 
point  is  observed  .  Singular  behavior  with  very  inte¬ 
resting  results  including  a  sudden  shift  of  breakdown 
point  is  observed  at  300  Hz  .  Finally  a  more  effective 
relaminar izat ion  than  observed  at  lower  Reynolds  num¬ 
ber  occurs  above  300  Hz  leading  to  a  completely  stabi¬ 
lized  laminar  flow  at  frequencies  above  380  Hz  . 

c.  In  the  case  of  turbulent  flow  ,  (Re=2640)  ,  no  notice¬ 
able  effect  is  observed  up  to  a  frequency  of  100  Hz. 
Above  this  frequency  the  breakdown  point  begins  to 
move  gradually  downstream  up  to  250  Hz  ,  returns  back 


upstream  to  a  frequency  of  333  Hz  ,  and  becomes  almost 
completely  unaffected  by  the  pressure  perturbations 
at  higher  frequencies. 

d.  Changing  the  waveform  of  the  disturbance  produced 
only  a  slight  deviation  in  the  average  breakdown 
point,  however  a  markedly  different  patern  of  droplet 
formation  was  observed. 

e.  Although  it  was  not  possible  to  record  the  phenomena 
photographically  for  inclusion  here  ,  it  was  possible 
using  the  phase  angle  control  incorporated  in  the 
stroboscope  to  observe  the  entire  development  of  the 
droplet  formation  visually.  This  process  was  found  to 
be  a  complex  sequence  of  breakup  and  re-agglomeration. 


Figure  34.  Photos  Illustrating  Jet  Breakdown  and  Oscilloscope 
Trace  of  Pressure  at  Re=2640  and  f=250  Hz 
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APPENDIX 


CALCULATIONS  AND  DATA  REDUCTION 


A.  CALCULATION  OF  CRITICAL  REYNOLDS  NUMBER 

A  visual  observation  of  the  flow  was  used  in  order  to 
determine  the  critical  Reynolds  number  .  Adjusting  the  flow 
by  the  means  of  the  needle  valve  ,  installed  at  the  entrance 
of  the  plenum  chamber  ,  a  laminar  jet  was  observed  up  to  the 
flow  meter  indication  itk  =  25  .  Above  this  value  a  turbulent 
flow  resulted. 

Using  the  equation  of  continuity  : 

m  =  pAU  or  U  =  m  /  pA  (Eq.  A-l) 

where  s 

m  =  the  fluid  mass  rate. 

U  =  the  mean  flow  velocity. 

p  =  62.34  lbm  /ft3  ,  the  water  density  (at  60°  F). 

A  =  the  cross  section  area  of  the  pipe. 

Also: 

A  =  Ttd2  14  =  (t r/  4 )  [( 3  /  16 )  ( 1  /  12 )]  2  (ft2) 
or:  A  =  1.9175  .10-4  (ft2) 
and : 

m  =  (.6(25/100) J(gal/min)  or  m  =  .15  gal/min 
or:  m  =  ( . 15  )  ( 500 . 8 )  ( 1/ 3600  )  (lbm/sec) 

or:  m  =  .020867  (lbm/sec) 

Substituting  the  above  values  into  Eq.  A-l  we  get: 


020867 


U  - - (ft/sec) 

(62.34Hl.9175  10-4 

or:  U  =  1.74562  (ft/sec) 

The  Reynolds  number  is  given  by: 

Re  =  (pUd)  /  (pgc)  (Eq.  A-2) 

where: 

p  =  2.4  .10"5  (lbf-sec/f t2) 

and : 

gc  =  32.2  ( 1 bm-f t )/(lbf-sec2) 

Substituting  the  above  values  into  (Eq.  A-2)  we  get: 

(62.34) (1 . 74562)(l/64) 

(Re)  =  - - - 

(2.4  10_5)(32.2) 

or:  (Re)  =  2200 
cr 

The  value  obtained  may  be  considered  reasonable  taking 
into  account  that  is  the  most  common  value  found  in  the  Hy¬ 
drodynamics  Bibliography.  The  numerical  values  used  in 
the  above  calculations  were  obtained  from  Ref.  25  . 

B.  CALCULATION  OF  INDICATED  PRESSURE  IN  LAMINAR  FLOW 

Refering  to  Fig.  B-l,  where  the  exact  dimensions  used 
are  indicated,  we  apply  the  Energy  Equation: 

P1  U12  e  p2  u2  2  8 

—  +  b.— +  —  z.  =  — -  +  b9— =—  +  -  z9  +  h.  (Eq.  B-l) 

P  2SC  Sc  P  2gc  gc 
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Figure  B-l.  Sketch  of  the  Apparatus. 


where  j 

s  the  pressure  at  point  1. 

P2  s  the  pressure  at  point  2  (  P2  =  patm  =  14.7  psia) 

s  the  mean  velocity  of  the  fluid  at  point  1. 

U2  :  the  mean  velocity  of  the  fluid  at  point  2. 

=  t>2  =  2  s  the  energy  correction  factor  for  laminar 

flow  at  points  1  and  2  respectively. 

:  the  height  difference  from  the  reference  plane 
:  the  flow  losses  given  by  the  relation: 

1  U2 

h,  =  f -  (Eq.  B-2 ) 

1  d  2gc 
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where: 

f  =  64/Re  :  friction  factor  for  laminar  flow. 

1  :  the  total  pipe  length. 

Inserting  the  numerical  values  and  the  appropriate  di¬ 
mensions  from  Fig.  B-l  into  Eq.  B-2  we  get: 

hL  =  (64/2200) (94.416/ .1875)  [(1. 74562)^/2(32. 2)j 
or:  hj  =  .69335  (lbf-ft/lbm) 

Using  the  Continuity  Equation  refering  to  Detail  A  of 
Fig.  B-l  and  neglecting  the  minor  turning  losses  we  get: 

m  =  pAUx  =  pAU2  (Eq.  B-3) 

or:  U2/U1  =  Ax/A2  =  d^/dj2  =  <  .25)2/(  .1875)2  =  1.7777 

Therefore: 

U1  =  V1’777?  =  1  •  74562/1 , 7777  (ft/sec) 

or:  Uj  =  .982  (ft/see) 

From  the  Energy  Equation,  (Eq.  B-l)  we  may  write: 

—  -  P2>  =-  (n22-U!2)  +  +  8chl 

P  2 

Solving  for  AP  =  P1~P2  we  get! 

=  -  [  -  (OjS-U^)  ♦  gU2-Zl>  ♦  (  B_4) 

gc  2 

Substituting  numerical  values  into  Eq.  B-4  we  get: 

4P  _  62 2_  ((lt74562)2.(,982)2]+32. 2i|+32. 2(. 69335)) 
32.201  2  12 


or : 


AP  =  125.18  ( lbf/ft2 ) 


or:  AP  =  .8693  (psig) 

or:  AP  =  1.775  (in.Hg)  at  60°  F 

The  latter  is  a  reasonable  result  compared  with  the 
mercury  manometer  which  indicated  pressure  Pjn<j  =  1.8  in.Hg 
suggesting  an  error  of  1.39  %  . 

The  curvature  of  the  pipe  does  not  create  a  large  effect 
in  friction  factor  for  this  Reynolds  number  which  is  in 
good  agreement  with  the  experimental  work  of  White  [Ref. 26]. 

C.  DYNAMIC  RESPONSE  OF  INDICATED  PRESSURE 

Since  the  pressure  perturbations  were  introduced  in 
the  upstream  plenum  chamber  and  measured  with  a  pressure 
transducer  through  a  pressure  tap,  we  may  reasonably  inquire 
into  the  accuracy  of  the  pressure  traces  measured.  In  order 
to  answer  this  question,  we  consider  the  dynamic  response  of 
the  system  under  the  assumption  that  the  fluid  pressure  is 
applied  uniformly  at  the  transducer  front  surface  . 

The  transient  response  of  a  pressure  measuring  instru¬ 
ment  is  dependent  on  two  factors: 

a.  The  response  of  the  transducer  element  that  senses 
ses  the  applied  pressure. 

b.  The  response  of  the  pressure  transmitting  fluid 
contained  in  the  connecting  tubing. 


This  latter  factor  is  frequently  the  one  which  deter¬ 
mines  the  overall  frequency  response  of  a  pressure  measure¬ 
ment  system  ,  and  eventually  calibration  must  be  relied  upon 
for  determining  this  response. 

An  estimate  of  this  behavior  may  be  obtained  from  the 
following  analysis  [Ref.  27]. 

Consider  the  system  shown  in  Fig.  C-l.  The  fluctuating 
pressure  has  a  frequency,  ,  and  an  amplitude,  PQ,  and  is 
impressed  on  a  tube  of  length  L,  and  radius  r.  At  the  end  of 
this  tube  is  a  chamber  of  volume  V,  where  connection  with 
the  pressure  transducer  is  made  .  We  seek  an  expression  for 
the  pressure  indicated  by  the  transducer,  PQ. 


Figure  C-l.  Schematic  Diagram  of  Pressure 
Transducer  Connection. 


The  mass  of  the  fluid  oscillates  under  the  influence 
of  fluid  friction  in  the  tube  which  tends  to  dampen  the  mo¬ 
tion.  If  the  conventional  formula  for  laminar  friction  resi¬ 
stance  in  tube  flow  is  used  to  represent  this  friction  ,  the 
resulting  expression  for  the  pressure  amplitude  ratio  is: 


[■£] =  {(l  ‘  <-^’)2  +  4h2^>2}-1/2  «*•  C-D 


where ; 


(*)n  is  the  natural  frequency  given  by  the  expressions 


r3  Trr  2c2 


%  = 


4LV 


(Eq.  C-2) 


and  h  is  the  damping  ratio  given  by  the  expression: 


h  = 


2Hg„  .  f3LV 


per *  1/  tt 


(Eq.  03) 


In  the  above  relations,  Eqs.  C-2  and  C-3  : 
c  =  the  velocity  of  sound  in  the  fluid, 
p  =  the  dynamic  viscosity  of  the  fluid, 
p  =  the  fluid  density. 

The  phase  angle  for  the  pressure  signal  is  given  by: 

(Eq.  04) 


.  ,  -1  -  2h<-£) 

4>  s  tan  1  - _1 


1  -  (-^-)2 
1  v  u*' 


Introducing  the  physical  dimensions  of  the  present  appa¬ 
ratus,  refering  to  Fig.  C-l  we  have: 

L  =  3.5  in  L  =  0 .5  in 

r  =  0.25  in  r' =  0.25  in 

V  =  rtr 2L  =  rr(  .25/12)2(3.5/12)  =  6.8177  10~5  (ft3) 
c  =  (K/p)1/2 

or:  c  =  (144  (311  . 103  )  / 1 . 938]  1/2  =  4800  (ft/sec) 

p  =  62.34  lbm/ft3 

p  =  2.4  .10-5  lbf-sec/f t3  (at  60°F 
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Substituting  the  above  values  into  Eq.  O 2  we  get: 


3tt(. 25/12)2(4800)2 
Wn  4(3.5/12)  (6.8177  10-4) 


or:  to Q  =  10885  see 


-1 


Substituting  into  Eq.  03  we  get: 


h  = 


2(2.4  10”5  )  (32.2)  \/ 3(3.5/12)(6.8177  10‘4) 


(62.34) (4800)(  .25/12) 

.-6 


rr 


or:  h  =  7.87  10 

From  these  results  we  see  that  the  difference  between 
pressure  transducer  indication  and  actual  plenum  chamber 
pressure  is  negligible  for  the  working  frequency  range. 

For  example  at  u)  =  1000  Hz  ,  from  Eq.  Ol  we  get: 

2 


P 

P_ 


( 

- 

1000' 

2' 

1- 

l 

L 

10885 

. 

+  4(7.87  10-6 ) 2 


1000 


10885 


2-1/2 


or:  P/P^  =  1.0085 

o 

Therefore  the  pressure  indication  error  due  to  dynamic 
response  is  0.85  % 

We  can  get  for  the  phase  angle  of  the  pressure  signal, 
from  Eq.  04  : 

=  tan"1  -[2(7.87  lO-6 )(. 09187)  /  (1  -  .09187)] 


or:  4*  =  -8.355  10-5  rad 
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D.  CALIBRATION  DIAGRAMS 

In  order  to  compare  theoretical  and  experimental  re¬ 
sults  concerning  the  indicated  pressure  from  the  mercury  ma¬ 
nometer,  as  well  as  the  pressure  variations  from  the  digital 
multimeter,  it  was  necessary  to  calibrate  both  of  them. 

A  pressure  calibration  diagram  is  plotted  in  Fig.  D-l  where 
it  is  shown  to  be  linear. 

Since  the  Reynolds  number  is  a  function  of  flow  rate,  it 
was  easier  to  adjust  the  flow  more  accurately  with  the  aid 
of  a  pressure  indicating  digital  multimeter  than  the  rotame¬ 
ter.  Therefore  two  diagrams  were  created,  shown  in  Figs.D-2 
and  D-3  where  are  plotted  the  indicated  flow  rate  versus  the 
indicated  pressure  ,  expressed  in  inches  Hg  and  D.C  Volts 
respectively. 


9 


P(IN.HG)  VS  P(D.C  VOLT) 


Figure  D-l.  Pressure  Calibration  Diagram 


FLOW  RATE  VS  PflN.HG) 


1.7  9.7  9.7  4.7  6.7  9.7  7.7  6.7  9.7  10.7  11.7 


PRESSURE  INDICATED  (IN.HG) 

Figure  D-2.  Flow  Rate  vs.  Pressure  Indicated  (in.Hg)  Diagram 
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